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ABSTRACT 


Reflected  shock  waves  were  used  to  heat  NgO/CgHg/CO/Ar  and  NgO/CH^/CO/Ar 
mixtures  to  temperatures  of  1875-2855  K  at  total  concentrations  of  5  x  10  cm'  . 
Oxygen  atom  production  was  monitored  via  the  flame  band  emission  at  450  nm,  and 
C02  production  was  observed  at  4.27  pm.  These  data  were  then  compared  to  the  re¬ 
sults  of  numerical  Integration  of  the  rate  equations,  using  rate  constant  and 
mechanistic  information  obtained  in  earlier  studies  of  H2  and  CH20  oxidation 
under  similar  conditions.  For  the  CgHg  mixture,  it  was  possible  to  achieve  good 
agreement  between  these  calculations  and  the  observed  data  using  only  one  addi¬ 
tional  reaction 

0  +  CH3  -►  CH20  +  OH. 

The  calculations  were  insensitive  to  the  precise  value  used,  but  the  results  are 
consistent  with  recent  high  temperature  literature  values.  For  the  CH^  mixture, 
the  best  agreement  was  achieved  using  a  recent  high  temperature  value  for 

0  +  CH4  -  CH3  +  OH. 

The  data  indicated  that  this  reaction  has  a  markedly  non-Arrhenius  rate  constant; 
use  of  a  value  based  upon  low  temperature  data  gives  results  inconsistent  with 
those  observed.  The  CH^  system  was  Insensitive  to  any  other  methane  reaction 
rate  constant.  There  was  some  experimental  evidence  to  suggest  an  unsuspected 
complexity  in  methyl  radical  decay  channels  at  the  lowest  temperatures  observed. 

The  mechanism/rate  constant  combination  used  here  was  then  applied  to  the 
analysis  of  literature  data  for  CH3  oxidation  by  02.  This  analysis  was  designed 
to  obtain  a  value  for  the  rate  constant  of 

CH3  +  02  ♦  CH20  +  OH. 

Results  Indicated  that  the  most  recent  high  temperature  studies  yield  values 
consistent  with  the  present  analysis. 
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INTR0DUCT10N 

This  paper  presents  the  results  of  the  next  step  in  our  study  of  klnetlcally 
simplified  oxidation  systems  [1*2].  In  an  attempt  to  more  accurately  characterize 
certain  high  temperature  rate  constants,  we  have  studied  systems  where  oxygen  atoms 
were  produced  at  well  defined  rates  via  NgO  dissociation;  the  oxygen  atom  profile 
was  then  followed  as  the  atoms  reacted  with  added  fuel  molecules.  In  principle 
such  an  approach  should  yield  reliable  rate  constant  information  for  the  reaction 
between  oxygen  atoms  and  the  fuel  molecule  of  interest.  In  practice  such  an  as¬ 
signment  can  be  quite  difficult  due  to  the- inevitable  presence  of  competing  reac¬ 
tions  of  various  types.  To  attempt  to  circumvent  this  problem,  this  approach  was 
first  applied  to  the  well -characterized  Hg  system  [1].  Although  it  might  appear 
that  the  next  step  would  be  to  study  methane,  it  soon  became  apparent  that  addi¬ 
tional  information  was  required  on  the  reactions  of  formaldehyde.  The  recent  com¬ 
pletion  of  such  a  study  [2]  suggested  this  system  was  sufficiently  well  character¬ 
ized  to  proceed  to  a  more  complex  system. 

Prior  to  study  of  the  practically  important  CH^  system,  it  was  thought  prudent 
to  study  a  system  involving  methyl  radicals.  The  first  phase  of  the  present  study 
focused  upon  reactions  of  methyl  radicals  produced  via  very  rapid  dissociation  of 
ethane.  This  approach  is  analogous  to  that  used  earlier  [3,4]  to  study  the  CH3/O2 
system.  The  added  advantage  of  the  CgHg/NgO  system  is  that  one  need  not  account 
for  reactions  of  Og  and  thus  the  analysis  is  more  straightforward,  particularly 
in  light  of  the  earlier  work  on  Hg  and  CHgO  under  similar  conditions. 

The  results  of  this  analysis  were  then  In  hand  when  the  ^O/CH^  data  were 
analyzed.  Here  the  prime  purpose  was  to  see  if  the  results  obtained  were  con¬ 
sistent  with  the  recent  work  [5,6]  which  indicated  that  the  rate  constants  for 


the  reactions  of  0  and  OH  with  CH^  had  a  non-Arrhenius  temperature  dependence. 
Such  Information  Is  particularly  important  with  respect  to  the  proper  choice 
of  rate  constants  for  use  In  modeling  studies  of  combustion  processes. 

The  results  described  herein  suggest  that  It  is  possible  in  most  Instances 
to  obtain  very  good  agreement  between  calculated  and  observed  profiles  using 
the  "best"  literature  values  for  the  pertinent  rate  constants.  However  there 
Is  some  evidence  to  suggest  that  there  may  be  an  additional  methyl  radical  decay 
route  at  the  lowest  temperatures  studied.  These  results  were  then  extended  to  in 
dude  an  analysis  of  the  CH3/02  system  with  particular  emphasis  upon  assignment 
of  a  high  temperature  rate  constant  for  the  reaction  between  methyl  radicals  and 
molecular  oxygen.  This  analysis  indicated  that  the  best  choice  here  was  consis¬ 
tent  with  that  recently  obtained  from  studies  of  more  complex  systems. 


-4- 


EXPERIMENTAL 

The  7.6  cm  shock  tube,  gas  handling  system,  and  optical  configuration  have 
been  described  previously  [7].  Infrared  emissions  were  collimated  by  two  slits 
1.5  nm  wide  and  5  mm  high  placed  50  mm  apart.  After  passing  through  a 
filter,  the  emissions  were  focused  upon  a  liquid-nitrogen-cooled  indium  antimonlde 
detector;  the  detector-preamp  combination  had  a  nominal  rise  time  of  1  us.  C02 
was  monitored  with  a  4.27  ym  interference  filter  (FWHM  =  0.18  ym).  The  flame- 
band  emissions  were  monitored  by  a  RCA  1P28A/V1  photomultiplier  mounted  behind  a 
450  nm  Interference  filter  (FWHM  =  6.5  nm)and  two  slits  1  mm  wide  by  2  mm  high 
which  were  38  mm  apart.  The  measured  response  time  of  the  photomultiplier/pre¬ 
amp  system  was  less  than  2  ys.  Data  were  collected  with  either  a  Biomation  Model 
805  transient  recorder  or  the  system  previously  described  [8],  Biomation  data 
were  collected  at  0.5  ys  Intervals  while  the  other  device  collected  data  at  10.0 
ys  Intervals. 

Gases  used  were  AIRC0  Ar  (99.9995%)  and  N20  (99.995%),  Matheson  CO  (99.99%), 
C2Hg  (99.96%)  and  CH^  (99.99%).  The  CO  was  slowly  passed  through  a  coil  of  copper 
tubing  at  77 K  before  addition  to  the  vacuum  line.  Other  gases  were  used  as  supplied. 

The  shock  tube  was  pumped  down  to  ^  2  mPa  and  the  observed  leak-outgasslng 
rate  was  usually  near  2  mPa/mln.  The  tube  was  Isolated  from  the  pumps  for  approx¬ 
imately  one  minute  prior  to  shock  initiation;  the  nominal  background  pressure  was 
near  4  mPa.  Since  the  test  section  was  pressurized  to  2.7  kPa  with  the  mixtures 
of  interest,  the  background  Impurities  for  the  shock  tube  were  present  in  the 
low  ppm  range  In  the  shocked  mixtures.  Mylar  diaphragms  and  helium  driver  gas 
were  used  throughout. 
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Temperatures behind  the  reflected  shocks  were  computed  In  the  usual  way  from 
measured  Incident  shock  velocities.  Reflected  shock  pressures  were  measured  with 
a  fast  response  pressure  transducer  and  were  always  found  to  be  In  good  agreement 
with  those  calculated.  Furthermore,  the  reflected  shock  pressure  was  always  con¬ 
stant  over  the  time  interval  that  data  were  collected;  as  a  result,  no  attempt  was 
made  to  correct  for  non-ideal  effects. 
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RESULTS 

Flame-band  Signal 

Information  about  oxygen  atom  production  was  obtained  by  observing  the 
flame-band  emissions  at  450  nm  in  systems  containing  added  CO.  The  procedure 
Is  completely  analogous  to  that  used  previously  [1].  Mixtures  studied  are 
listed  In  Table  1.  As  In  the  earlier  work,  additional  mixtures  were  prepared 
which  were  similar  to  those  listed  except  that  CO  was  omitted.  Here  the  time- 
resolved  emissions  at  450  nm  were  recorded  and  used  as  background  for  the  ex¬ 
periments  of  interest.  The  calibration  factor  relating  the  flame-band  inten¬ 
sity  to  the  product  of  the  CO  and  0  concentrations  was  the  same  as  that  used 
earlier.  Additional  checks  of  this  factor  during  the  course  of  these  experi¬ 
ments  verified  the  value.  Data  were  collected  at  0.5  us  intervals.  For  con¬ 
venience  in  the  data  reduction  process,  times  were  measured  relative  to  tQ,  the  time  at 
which  the  reflected  shock  was  calculated  to  be  first  visible  at  the  detector.  This 
point  was  determined  by  simple  geometrical  optics  arguments  using  the  collimating 
slit  dimensions  employed  for  these  experiments.  Typically  tQ  is  ^  5  ys  earlier 
than  tmp,  the  time  of  shock  passage  of  the  window  midpoint. 

The  results  of  typical  experiments  are  shown  in  Fig.  1.  Here  three  para¬ 
meters  were  used  to  characterize  the  signal.  The  first,  t^2»  designates  the 
time,  relative  to  tQ,  when  the  signal  reaches  one-half  its  maximum  value.  Simi¬ 
larly,  is  the  time  when  the  signal  reaches  three-fourths  the  maximum.  Vm 
specifies  the  maximum  signal  observed.  Three  factors  were  explicitly  considered 
In  an  error  analysis  of  these  data:  (1)  Measurement  errors,  including  noise,  cal¬ 
ibration  errors,  and  errors  in  background  corrections,  (2)  errors  of  ±2  ys  in 
location  of  tQ,  and  (3)  errors  resulting  from  temperature  uncertainties.  More 
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explicit  details  of  the  procedures  used  In  (2)  and  (3)  can  be  found  In  Ref.  2. 
These  errors  were  combined  In  the  standard  way  [9]  to  estimate  the  overall  un¬ 
certainty.  These  parameters  are  shown  in  Figs.  2  and  3  and  listed  In  Table  1. 

Note  that  the  scatter  observed  appears  to  be  consistent  with  the  estimated  error 
bars.  These  parameters  appear  to  adequately  characterize  the  formation  of  oxygen 
atoms  In  these  systems.  Although  the  t^g  values  are  comparable  In  Mixtures  A 
and  B,  there  seem  to  be  small  but  real  differences  in  t and  Vm;  the  oxygen 
atoms  appear  to  Increase  somewhat  more  slowly  in  the  CgHg  mixture  and  ultimately 
achieve  a  higher  steady  state  concentration. 

COg  Production 

Infrared  emission  at  4.27  pm  was  observed  at  10.0  ys  Intervals.  Here  times 
were  measured  relative  to  the  passage  of  the  window  midpoint  by  the  reflected 
shock  (tmp).  Use  of  this  procedure  simplified  comparisons  of  calculations  and 
observations.  The  slower  sampling  rate  here  minimized  perturbations  due  to  the 
finite  window  width. 

The  observed  signal  was  corrected  for  emissions  from  NgO  and  CO  using  pro¬ 
cedures  described  previously  [2].  Typical  data  are  shown  In  Fig.  4,  Three  para¬ 
meters  were  used  to  characterize  these  profiles:  tQ  5  Indicates  the  time  (rela- 
tlve  to  tmp)  that  [COg]  =  5x10  /cm  ;  this  was  considered  to  be  the  lowest  con¬ 
centration  at  which  reliable  data  could  be  obtained.  t|  5  and  tg  g  Indicate  the 
times  [COg]  ■  1.5  x  1016/cm3  and  2.5  x  lO^/cm3,  respectively.  The  error  analysis 
procedure  used  here  was  similar  to  that  used  for  the  flame-band  data.  In  some 
cases  very  large  error  bars  result  from  the  fact  that,  after  the  "knee"  In  the 
profile,  small  concentration  errors  result  In  large  time  errors.  Plots  of  these 
parameters  are  shown  in  Figs,  5  and  6;  they  are  listed  in  Table  1,  Again  it  Is 
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seen  that  the  observed  scatter  is  consistent  with  the  calculated  error  bars. 

Note  that  the  ($2  parameters  for  both  mixtures  are  similar.  The  minimum  seen 
in  the  Figs.  5(c)  and  6(c)  results  from  diminished  CO2  production  at  the  highest 
temperatures.  For  example,  the  concentration  of  C02  observed  at  500  us  decreases 

from  4  x  1016/cm3  near  2000  K  to  less  than  3  x  10^6/cm3  at  2800  K. 

i 
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DISCUSSION 

As  outlined  In  the  Introduction,  the  thrust  of  this  program  has  been  to  at¬ 
tempt  to  analyze  relatively  simple  kinetic  systems  of  importance  in  combustion. 

The  emphasis  has  been  to  minimize  the. kinetic  ambiguity  by  proceeding  very  slowly 
toward  systems  of  increasing  complexity.  In  this  manner,  the  number  of  potential 
variables  is  reduced  to  the  point  where  a  meaningful  analysis  should  be  possible. 

In  practice,  this  has  been  achieved  by  utilizing  all  of  the  information  from  past 
work  as  fixed  quantities  not  subject  to  modification  in  an  attempt  to  fit  the  new 
data.  Of  course,  some  modifications  may  ultimately  be  required,  e.g.,  the  change  - 
in  the  rate  constant  assignment  for  H  +  NgO  a  N2  +  OH  [1,2],  but  this  approach 
serves  to  provide  a  much  more  constrained  system  in  which  fortuitous  agreement  of 
observations  and  calculations  would  appear  to  be  more  unlikely. 

The  likely  mechanism  for  the  N20/C2Hg/C0/Ar  system  (Mixture  A)  (assuming  suf¬ 
ficiently  rapid  C2Hg  dissociation)  is  listed  in  Table  2  as  Reactions  1-19.  The 
most  important  feature  of  this  mechanism  is  that,  in  spite  of  its  apparent  com¬ 
plexity,  information  is  available  about  eighteen  of  these  reactions  from  studies 
of  simpler  systems  with  the  same  experimental  apparatus  used  in  this  work.  Given 
the  fact  that  good  agreement  between  experiments  and  calculations  were  achieved 
there,  one  can  consider  the  present  analysis  to  consist  of  simply  analyzing  the 
effect  of  the  one  additional  reaction,  namely  that  between  oxygen  atoms  and 
methyl  radicals  (Reaction  5). 

In  the  calculations  with  this  mechanism,  provisions  were  made  for  both  the 
finite  observation  window  width  and  detector  response  times  as  in  the  earlier 
work  [2].  The  flame-band  calculations  included  the  calibration  factor  so  that  the 
computed  output  was  in  effect  an  absolute  oxygen  atom  concentrati on-time  profile; 


no  normalization  factors  were  used.  The  slit  corrected  flame-band  output  could 
then  be  directly  compared  to  the  observations  with  t  as  the  reference  time.  For 
C02  comparisons,  slit  effects  were  ignored  since  data  were  obtained  at  10.0  us 
intervals.  Here  t  was  the  reference  time, and  again  it  was  possible  to  directly 
compare  calculated  and  observed  profiles  using  the  experimental  C02  calibration 
experiments  [1,2]. 

Comparisons  were  done  at  selected  temperatures.  In  most  cases  the  experimental 
parameters  were  taken  directly  from  Table  1;  in  a  few  instances  the  values  were  ad¬ 
justed  slightly  to  better  conform  to  other  data  in  that  temperature  range.  In  cases 
where  adjustments  were  made,  the  magnitude  of  the  shift  was  always  smaller  than  the 
estimated  error  bar.  These  representative  experimental  parameters  for  Mixture  A, 
with  the  estimated  errors,  are  listed  in  Table  3. 

The  calculated  profiles  were  treated  identically  to  the  experimental  ones  to 
obtain  calculated  values  of  the  reaction  par?ieters.  The  mechanism  and  rate  con¬ 
stants  listed  in  Table  2  [Reactions  (1)-(19)]  were  designated  Set  1  and  the  results 
are  listed  in  Table  3.  These  results  are  also  illustrated  in  Figs.  2  and  5,  and 
the  total  profile  is  shown  in  Fig.  1(a)  and  4(a).  The  agreement  seen  here  is  quite 
reasonable.  The  possible  problem  areas  include  calculated  ty2  values  which  are 
somewhat  long,  a  calculated  temperature  dependence  of  Vm  somewhat  higher  than  ob¬ 
served,  and  a  tendency  for  C02  production  to  be  a  little  slow  at  long  times.  How¬ 
ever,  these  problems  are  quite  minor  in  light  of  the  diverse  data,  and  the  fact 
that  the  comparisons  are  made  on  an  absolute  concentration-time  basis.  Particularly 
encouraging  was  the  fact  that  no  adjustments  of  any  rate  constants  assigned  in  the 
earlier  work  [1,2]  were  needed  to  achieve  this  fit. 

The  justification  for  the  assumption  that  C2Hg  dissociation  is  sufficiently 
rapid  so  that  its  precise  rate  was  unimportant  was  checked  in  a  series  of  calculations 
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using  20  reactions.  Here  C2Hg  dissociation  was  included  using  a  rate  constant 
equal  to  one-half  that  recently  reported  by  Olson,  et  al.  [22].  (The  lower  value 
was  used  to  provide  a  more  rigorous  test  of  the  assumption.)  This  calculation  pro¬ 
duced  results  virtually  identical  to  those  obtained  in  Set  1. 

The  sensitivity  of  the  results  to  variations  in  kg  were  checked  in  Sets  2  and 
3  where  the  rate  constant  was  varied  both  directions  by  a  factor  of  two.  Note 
the  results  are  surprisingly  insensitive  to  this  variation.  These  experiments 
are  not  capable  of  providing  a  more  reliable  high  temperature  estimate  of  this 
rate  constant.  It  would  appear  that  the  best  value  would  be  in  the  range  kg  = 

1-2  x  10’10  cm3  s’1.  This  value  is  consistent  with  those  frequently  used  at 
high  temperatures  [11,14]  as  well  as  the  recent  low  temperature  measurements 
[23,24]. 

Attempts  were  then  made  to  see  if  the  remaining  disagreements  between  cal¬ 
culations  and  observations  could  be  related  to  other  reactions  of  CH^.  This  was 
checked  by  additional  calculations  using  a  23  reaction  mechanism  in  which  the 
following  reactions  were  included: 


CH3  +  M 

=  CH2  +  H  +  M 

(24) 

ch3  +  o2 

=  CH20  +  OH 

(25) 

CH3  +  H 

=  ch2  +  h2 

(26) 

CH3  +  OH 

=  ch2o  +  h2 

(27) 

The  rate  constants  for  these  reactions  were  taken  from  Olson  and  Gardiner  [25]. 
The  results  are  listed  as  Set  4  in  Table  3.  Again,  there  are  only  marginal  dif¬ 
ferences  in  the  calculated  parameters,  suggesting  these  reactions  are  of  little 
consequence  under  these  conditions.  The  slightly  slower  C02  production  seen  here 
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1s  probably  due  to  the  fact  that  some  of  the  CH3  Is  converted  to  CH2  which  Is 
"unreactive11  In  this  mechanism  and  serves  as  a  sink  for  hydrogen  atoms. 

In  spite  of  the  Inability  of  the  present  analysis  to  more  precisely  deter¬ 
mine  kg,  the  overall  agreement  here  presents  encouraging  evidence  that  the  earlier 
analyses  of  the  H2  and  CH20  systems  [1,2]  were  reasonable.  In  retrospect,  this 
system  simply  provides  a  method  for  rapid  generation  of  CH20  via  Reaction  (5). 

This  alternative  method  of  generating  CH20  suggests  that  the  formaldehyde  con¬ 
centration  in  the  earlier  work  was  indeed  in  the  expected  ranges  and  there  were 
no  significant  adsorption  problems  there. 

Calculations  with  Mixture  B  were  initiated  using  all  of  the  reactions  shown 
in  Table  2.  It  is  important  to  note  that  the  four  reactions  involving  CH^  all 
have  rate  constants  which  reflect  the  recent  high  temperature  data  available  on 
these  systems.  The  three  radical  reactions  with  methane  have  rate  constants  well 
above  a  simple  Arrhenius  extrapolation  of  low  temperature  data.  The  results  with 
this  mechanism/rate  constant  combination  are  shown  in  Table  4  (Set  I)  as  well  as 
in  Figs.  3  and  6.  A  representative  total  profile  calculation  is  shown  in  Fig.  1(b) 
and  4(b).  In  most  cases,  it  can  be  seen  that  there  is  excellent  agreement  between 
calculated  and  observed  results.  The  one  real  problem  concerns  at  2085K.  Here 
the  calculation  shows  a  distinct  "hitch"  just  after  t^2  where  there  is  a  dramatic 
decrease  in  the  rate  of  0-atom  production  for  a  short  time.  There  is  no  evidence 
for  such  behavior  in  the  observed  profiles.  However,  it  must  be  remembered  that 
the  flame-band  data  were  corrected  for  background  emissions  and  it  is  possible,  but 
unlikely,  that  the  discrepancy  reflects  an  artifact  produced  during  the  data  reduc¬ 
tion  process. 

Calculations  in  which  k21,  k22,  and  k23  were  varied  by  factors  of  two  each 
direction  indicated  the  profiles  are  not  sensitive  to  these  rate  constants.  The 
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effect  of  varying  k2Q  Is  shown  in  Sets  II  and  III.  Increases  in  k2Q  (Set  II) 
cause  obvious  difficulties  with  the  t^2  flame-band  data  while  slightly  Improving 
the  C02  results.  Decreases  in  k2Q  (Set  III)  have  a  mixed  effect  on  the  flame- 
band  data  while  causing  C02  production  to  be  slower  than  observed.  It  would  ap¬ 
pear  that  the  Roth-Just  value  used  in  Set  I  is  probably  a  reasonable  compromise 
here.  Set  IV  illustrates  the  effect  of  using  a  k2Q  value  based  upon  low  tempera¬ 
ture  flow  system  data  [26].  Although  the  ty2  and  t3/,4  values  at  2085K  appear  to 
be  a  distinct  improvement,  the  complete  profile  exhibits  a  significant  dip  just 
beyond  t^,^  which  is  incompatible  with  the  experimental  observations.  Further¬ 
more  note  that  there  is  uniform  disagreement  here  with  the  C02  data.  It  would 
appear  the  present  data  require  a  k2Q  value  well  above  what  one  would  expect  from 
the  low  temperature  flow  work. 

The  tj/4  discrepancy  at  2085K  suggests  there  may  be  some  defects  in  the 
mechanism/ rate  constant  combination  used  here.  Detailed  analysis  of  the  low 
temperature  results  on  Mixtures  A  and  B  revealed  that  both  calculations  showed 
the  hitch  in  the  flame-band  production.  The  usual  integration  precautions  (i.e. , 
decreasing  step  size,  etc.)  indicated  the  problem  was  not  a  computational  artifact. 
Instead  it  resulted  from  the  peaking  of  the  oxygen  atom  decay  rate  due  to  the  0 
+  CH^  reaction.  This  phenomenon,  coupled  with  an  essentially  constant  oxygen 
atom  production  rate,  causes  the  net  production  rate  to  pass  through  a  distinct 
minima.  The  problem  was  not  obvious  in  analysis  of  Mixture  A  since  the  hitch 
occurred  at  times  well  before  ty2  there.  However,  the  fact  that  CH3  is  produced 
at  later  times  in  Mixture  B  cause  this  hitch  to  be  more  obvious.  (Note  in  Set  IV 
that  a  lower  k2Q  value  further  delays  and  exacerbates  this  feature,) 

Assuming  that  the  flame-band  data  is  correct,  this  discrepancy  between  ob- 
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servatlons  and  calculations  suggests  that  there  may  be  an  alternate  CH^  decay 
route.  Such  a  route  would  decrease  0-atom  decay  via  0  +  CHj  and  thus  tend  to 
minimize  the  hitch.  Additional  methyl  reactions  were  included  in  another  cal¬ 
culation  (Set  V)  to  see  if  any  effects  were  observed.  Here  Reactions  ( 24 ) - ( 27 ) 
were  included  with  the  rate  constants  quoted  by  Olson  and  Gardiner  [25].  Note 
In  Table  4  that  these  adjustments  did  not  improve  the  calculated  parameters. 

It  would  appear  that  some  other  methyl  radical  reaction  must  be  responsible, 
but  it  is  not  clear  at  this  time  what  reaction(s)  could  account  for  this  dis¬ 
crepancy.  Note  that  the  presence  of  an  alternate  decay  mechanism  for  methyl  rad¬ 
icals  would  also  improve  the  agreement  in  the  t y2  measurements  in  Mixture  A. 

Here  suppression  of  the  oxygen  atom  "hitch"  would  result  in  a  lower  ty2  value. 

In  better  agreement  with  the  observations. 

It  is  Important  to  recognize  that  the  present  indications  of  some  unexpected 
complexities  in  the  mechanism  of  methyl  radical  decay  are  based  on  fragmentary 
evidence  collected  at  one  extreme  of  the  temperature  range  observed.  It  is  clear 
that  further  work  is  required  to  resolve  the  problem  posed  here.  An  obvious  route 
for  future  experimentation  would  involve  direct  observation  of  methyl  radicals  utiliz 
Ing  the  high  temperature  spectroscopic  data  recently  obtained  by  Glanzer,et  al.,[27]. 
It  would  appear  that  the  best  information  could  be  obtained  at  the  lower  tempera¬ 
tures  where  there  is  a  much  better  balance  between  oxygen  atom  production  and  de¬ 
cay  channels.  Given  the  high  activation  energy  for  NgO  dissociation,  the  pro¬ 
duction  route  Is  sufficiently  large  at  higher  temperatures  to  mask  any  subtle 
changes  In  the  decay  mechanism. 

In  spite  of  the  potential  difficulty  concerning  the  low  temperature  methyl 
decay,  it  is  apparent  that  the  mechanism/rate  constant  combination  given  in  Table 


■J.  '  W  '  — '  JL*  "IT*  — ’ 


-15- 


2  can  accurately  describe  most  of  the  data.  In  particular.  It  Is  apparent  that 
one  must  assume  a  non-Arrhenius  temperature  dependence  for  the  reaction  between 
oxygen  atoms  and  methane  to  fit  the  present  data.  This  conclusion  is  not  con¬ 
tingent  upon  the  other  methane  rate  constant  assignments  since  the  calculations 
were  insensitive  to  these  values. 

It  is  possible  to  use  the  mechanism/rate  constant  combination  developed  in 
the  earlier  work  [1,2]  and  confirmed  by  the  agreement  achieved  in  the  present 
study  to  analyze  the  data  of  Izod,  Kistiakowsky,  and  Matsuda  [3]  with  particular 
emphasis  upon  determination  of  the  rate  constant  for  Reaction  (25),  a  most  con¬ 
troversial  reaction  [28,29].  I KM  used  azomethane  as  a  precursor  of  methyl  radi¬ 
cals  and  then  studied  the  oxidation  in  a  C0/02  environment  quite  similar  to  that 
used  in  our  CHgO  oxidation  study  [2].  Thus,  analysis  of  the  IKM  data  can  be  con¬ 
sidered  as  a  straightforward  extension  of  that  work  with  only  two  additional  reac¬ 
tions  to  consider: 


0  +  CH3 

-  CH20  +  H 

(  5) 

02  +  CH3 

-►  CH20  +  OH 

(25) 

Attention  was  limited  to  consideration  of  the  azomethane  data;  the  analysis  of 
the  ethane  data  in  [3]  would  be  dependent  upon  the  ethane  dissociation  rate  due 
to  the  lower  concentration  and  temperatures  in  that  work. 

In  addition  to  the  two  reactions  listed  above  the  calculations  included 
Reactions  (6)-(19)  of  Table  2  and 

CO  +  02  -  C02  +  0.  (28) 

The  results  listed  in  Set  A  used  the  rate  constants  in  Table  2,  k2g  from  Olson 
and  Gardiner  [25],  and  k2g  from  our  earlier  work  [2].  Although  IKM  also  reported 
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exponential  growth  rate  constants,  these  are  not  listed  In  Table  5  since  semi- 
log  plots  of  the  calculated  profiles  were  sufficiently  non-linear  to  make  as¬ 
signments  of  A  ambiguous.  Although  data  and  calculations  are  listed  at  1500K 
for  completeness,  no  serious  attempts  were  made  to  Interpret  this  Information 
since  so  many  of  the  rate  constants  used  were  determined  at  higher  temperatures, 
there  were  simply  too  many  potential  variables  to  permit  a  meaningful  analysis 
here. 

Of  particular  interest  in  the  Set  A  comparisons  is  the  good  agreement  in¬ 
dicated  for  t..  at  both  1750K  and  2150K.  Although  the  fit  for  a  was  not  as  good. 

It  appears  to  be  reasonable. 

The  sensitivity  of  these  calculations  was  then  checked  in  a  variety  of  ways. 
The  effect  of  methyl  radical  recombination  was  analyzed  by  using  the  recombination 
rate  used  by  IKM  (k  *  1.4  x  10"11  cm"3  $"^).  These  results  (Set  B)  showed  little 
effect  above  1500K;  use  of  a  more  recent  value  [27]  which  is  significantly  lower 
than  that  used  by  I KM  Indicated,  as  expected,  even  less  of  an  effect.  Variation 
of  kg  (0  +  CH^)  showed  virtually  no  effect,  but  substantial  changes  were  seen  when 
^25  (®2  +  ^*3)  was  var^e<*'  These  changes  are  illustrated  by  Sets  B  and  C  in  Table 
5.  Note  tj  is  quite  sensitive  to  k25;  a  decrease  of  a  factor  of  two  yields  values 
significantly  different  from  those  observed.  Values  somewhat  larger  than  that 
used  by  Gardiner  might  be  acceptable;  assignment  of  limits  is  difficult  here  since 
errors  are  not  assigned  in  the  published  data.  However,  it  does  appear  that  a 
value  similar  to  either  Olson  &  Gardiner's  [25]  or  Tabayashi  and  Bauer's  [28] 
would  be  appropriate.  It  is  significant  that  this  value  is  well  above  the  upper 
limit  reported  by  Baldwin  and  Golden  [29]  at  1220K. 

The  agreement  seen  in  Table  5  is  not  as  impressive  as  that  reported  by  IKM  in 
their  calculations.  However,  the  present  fit  was  achieved  with  fewer  adjustable 
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parameters  and  a  more  complete  mechanism.  It  is  interesting  to  note,  however, 
that  the  value  reported  by  IKM  for  kgj  is  quite  close  to  that  used  here.  At 
2000K,  the  present  value  is  only  40%  larger.  The  major  difference  in  the  mechan¬ 
ism  used  in  this  analysis  is  the  use  of  all  the  formaldehyde  decay  channels  rather 
than  just  the  oxygen  atom  reaction  used  by  IMK.  As  a  result  the  rate  constant 
for  this  reaction  is  appreciably  slower  in  this  work. 

i 

Of  course,  the  disclaimer  must  be  added  that  the  kgg  assignment  is  contin¬ 
gent  upon  the  proper  mechanistic  assignment.  If  there  is  an  alternate  methyl 
decay  route,  as  might  be  indicated  by  the  low  temperature  data  in  Mixtures  A  and 
B,  it  is  conceivable  that  this  rate  constant  assignment  could  be  too  high.  How¬ 
ever,  the  alternate  route  in  this  case  must  yield  COg  at  a  rate  comparable  to  that 
proposed  here  to  be  consistent  with  the  data. 

In  spite  of  the  remaining  ambiguity,  it  does  appear  that  the  present  method 
of  assigning  kgg  is  more  direct  than  previous  attempts  since  the  number  of  vari¬ 
ables  was  kept  to  a  minimum.  However,  it  is  obvious  that  additional  experiments 
with  otner  diagnostic  tools  are  required  to  remove  the  remaining  uncertainties  In 
even  these  "kinetically  simplified"  systems. 


CONCLUSIONS 


It  has  been  possible  to  accurately  describe  oxygen  atom  and  C02  production 
for  the  oxidation  of  CH^  and  CH^  by  NgO  over  an  extended  temperature  range.  It 
was  not  necessary  to  adjust  any  rate  constants  to  achieve  the  observed  fit.  The 
CH4  analysis  confirmed  the  non-Arrhenius  character  of  the  reaction  0  +  CH^  -*>  CH^ 

+  OH.  It  was  also  possible,  using  the  set  of  rate  constants  which  have  now  been 
shown  to  be  consistent  with  a  relatively  large  data  base,  to  analyze  some  earlier 
experiments  on  CHj  oxidation  by  02  to  obtain  an  estimate  of  the  rate  constant  for 
the  reaction  CH3  +  02  -*■  CH20  +  OH.  It  would  appear  that  this  reaction  is  important 
in  hydrocarbon  oxidation,  and  its  rate  constant  is  now  known  reasonably  well.  It 
is  noted  that  some  problem  areas  still  exist  even  in  these  kinetically  simplified 
systems,  and  more  work  is  required  here  before  one  can  confidently  proceed  to  the 
analysis  of  the  more  complicated  hydrocarbon  oxidation  systems. 
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TABLE  2.  MECHANISM  AND  RATE  CONSTANTS 


Rate  Constant3 


Reaction 

l0£in_A 

n 

E. 

Reference 

1 .  N20  +  M 

“  N2  ♦  0  +  M 

-9.57 

0 

n — 

216 

1 

2.  N20  +  0 

-  2N0 

-10.11 

0 

117 

10 

3.  N20  +  0 

*  N2  +  02 

-10.11 

0 

117 

10 

4.  N20  +  H 

*  N2  +  OH 

-8.82 

0 

92 

2 

5.  0  +  CH3 

»  CH20  +  H 

-10.00 

0 

0 

11 

6.  0  +  CH20 

*  HCO  +  OH 

-10.52 

0 

12.9 

12 

7.  H  +  CH20 

*  HCO  +  H2 

-9.26 

0 

44 

2 

8.  OH  +  CH20 

*  HCO  +  H20 

-10.90 

0 

0.6 

13 

9.  M  +  CH20 

-  HCO  +  H  +  M 

-7.26 

0 

339 

2 

10.  0  +  HCO 

=  CO  +  OH 

-9.77 

0 

0 

14 

11.  H  +  HCO 

*  CO  +  h2 

-9.48 

0 

0 

14 

12.  OH  +  HCO 

*  CO  +  H20 

-9.77 

0 

0 

14 

13.  M  +  HCO 

=  CO  +  H  +  M 

-9.59 

0 

61.4 

15 

14.  0  +  CO  +  M 

=  C02  +  M 

-33.80 

0 

0 

10 

15.  OH  +  CO 

=  C02  +  H 

-11.18 

0 

33.4 

16 

16.  0  +  H2 

»  OH  +  H 

-9.44 

0 

57.3 

17 

17.  OH  +  H2 

*  H20  +  H 

-10.44 

0 

21.5 

18 

18.  H  +  02 

=  OH  +  0 

-6.69 

-.91 

69.5 

19 

19.  0  +  H20 

*  20H 

-9.96 

0 

76.9 

20 

20.  0  +  CH4 

=  ch3+  OH 

-9.17 

0 

58.5 

5 

21 .  H  +  CH4 

=  ch3  +  h2 

-8,92 

0 

63.1 

21 

22.  OH  +  CH4 

=  Crt3  +  H20 

-17.59 

2.13 

10.2 

6 

23.  M  +  CH4 

=  CH3  +  H  +  M 

-6.11 

0 

390 

21 

Expressed  in  the  form  k  *  ATnexp(-E^/RT)  in  units  of  cm^,  molecule,  sec,  and  kJ. 


TABLE  3.  COMPARISON  OF  CALCULATIONS  AND  OBSERVATIONS  (MIXTURE  A) 


Parameter3 

I/K 

Observed 

Set  lb 

Set  2C 

Set  3d 

Set  4e 

Flame  band 

2020 

28±  5 

36 

36 

33 

34 

Vi  /o 

2410 

9±  2 

13 

14 

— 

14 

1/2 

2855 

7±  2 

9 

9 

9 

9 

2020 

47±13 

54 

52 

55 

53 

1  4 

2410 

17±  4 

19 

19 

— 

20 

3/4 

2855 

10±  3 

11 

11 

11 

11 

2020 

0.36±.06 

0.29 

0.29 

0.30 

0.30 

vm 

2410 

1 . 2  5±  .15 

1.28 

1.25 

-- 

1.32 

m 

2855 

2.83±.32 

3.29 

3.20 

3.38 

3.37 

C02 

2020 

23±  4 

26 

24 

28 

28 

t  , 

2410 

8±  2 

9 

8 

9 

.  0 

2855 

“ 

5 

4 

5 

5 

2020 

43±  5 

49 

47 

51 

52 

h.s 

2410 

15±  3 

18 

17 

19 

2855 

14±  3 

11 

11 

14 

13 

2020 

65±  6 

77 

77 

80 

82 

t-  c 

2410 

26±  4 

34 

33 

-- 

37 

2.5 

2855 

370±80 

509 

448 

560 

525 

aSee  TABLE  1  for  definitions. 

Reactions  1-19  (TABLE  2). 
cHere  kg  *  2x1  O'10. 
dHere  kg  *  5x1 O"11. 

eHere  Included  additional  CH3  reactions  (see  text). 


TABLE  4:  COMPARISON  OF  CALCULATIONS  AND  OBSERVATIONS  (MIXTURE  B) 


Parameter8 

I/K 

Observed 

Set  I1 

Flame  band 

2085 

17±  3 

18 

* 1/2 

2435 

10±  3 

12 

2840 

7±  2 

9 

2085 

25±  6 

43 

*3/4 

2435 

15±  5 

17 

2840 

81  3 

11 

2085 

0.351.05 

0.32 

2435 

1.201.13 

1.16 

2840 

2.581.26 

2.82 

co2 

1875 

601  4 

58 

f 

2085 

22i  4 

24 

2435 

81  2 

8 

2840 

- — 

5 

1875 

1151  7 

104 

t.  „ 

2085 

35i  6 

41 

*1.5 

2435 

14i  3 

16 

2840 

--- 

9 

1875 

1601  9 

152 

♦  „ 

2085 

56±  7 

61 

*2.5 

2435 

261  4 

28 

2840 

300153 

289 

Set  IIC 

Set  IIId 

Set  IVe 

Set  Vf 

27 

15 

15 

18 

mm 

mm 

9 

9 

9 

— 

43 

39 

21 

42 

“ 

-- 

mm 

mm 

11 

11 

11 

— 

0.30 

0.33 

0.37 

0.32 

-- 

-- 

— 

2.75 

2.89 

3.01 

•m 

MM 

21 

27 

33 

i  ♦ 
1  N 

~ 

•• 

mm 

4 

5 

5 

m  m 

38 

45 

51 

43 

-- 

mm 

-- 

9 

10 

10 

mm 

58 

65 

72 

Mm 

64 

-- 

MM 

MM 

46 

339 

410 

MM 

aSee  TABLE  1  for  definitions. 

Reactions  1  -  23  (TABLE  2) 

C*c20  *  1-36  x  10"9  exp(-58.5  kJ/RT). 

dk2Q  •  3.4  x  10’10  exp(-58.5  kJ/RT). 

•kgQ  »  3.5:;x  10’11  exp(-37.8  kJ/RT). 

Included  additional  CHj  reactions  (see  text). 


TABLE  5:  COMPARISON  OF  CALCULATIONS  AND  OBSERVATIONS 


(235  ppm 

(CH3)2N2, 

4.18%  CO,  2.13%  C 

>2  in  Argon) 

Parameter 

VK 

Observed8 

Set  Ab 

Set  Bc 

Set  Cd 

Set  De 

1500 

300 

180 

215 

128 

270 

1750 

no 

no 

115 

80 

160 

2150 

60 

64 

65 

45 

90 

1500 

2 

0.9 

0.9 

0.9 

0.9 

a9  1750 

3 

2.0 

2.0 

2.0 

2.0 

2150 

5 

4.1 

4.4 

4.0 

4.1 

aTaken  from  Figs.  2(a)  and  3(a)  of  Reference  3. 

^"Standard  Set";  see  text. 

cIncluded  CH3  +  CH3;  see  text. 

dk25  »  2.4  x  10"12  exp(-38  kJ/RT). 

ek25  »  6.0  x  10'13  exp (-38  kJ/RT). 

fT1me  [C02a  *1.6x  1013/cm3. 

^Constant  (late)  rate  of  C02  production. 


LEGENDS 


Fig.  1 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


Fig.  6 


Typical  Flame-band  Profiles,  (a)  Mixture  A,  T  *  2410K.  (b)  Mixture  B, 
T  -  2435K.  Open  circles  designate  every  fifth  data  point;  closed  cir¬ 
cles  are  the  calculated  results.  (Sets  1  and  I.) 


Flame-band  Parameters  for  Mixture  A.  (a)  Time  required  to  reach  one- 
half  maximum  voltage,  (b)  Time  required  to  reach  three-fourths  maxi¬ 
mum  voltage,  (c)  Maximum  voltage.  Open  symbols  are  experiments, 
closed  symbols  are  calculated  values.  (Set  1.)  All  times  are  relative 
to  tQ.  (See  text.) 


Flame-band  Parameters  for  Mixture  B.  Symbols  have  same  meaning  as 
In  Fig.  2.  Calculations  are  Set  I. 

Typical  C02  Profiles,  (a)  Mixture  A,  T  =  2410K.  (b)  Mixture  B, 

T  *  2435K.  Open  circles  are  observed;  closed  circles  are  calculated. 
(Sets  1  and  I.) 


16  3 

C02  Parameters  for  Mixture  A.  (a)  Time  for  [C02]  to  reach  0.5  x  10  /cm  . 
(b)  Time  to  reach  1.5  x  10^/cm  .  (c)  Time  to  reach  2.5  x  lO^/cm'*. 

Open  symbols  are  observed;  closed  symbols  are  calculated.  (Set  1.) 

All  times  measured  relative  to  tmp.  (See  text.) 


C02  Parameters  for  Mixture  B.  Symbols  have  same  meaning  as  In  Fig.  5. 
Calculations  are  Set  I. 
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additional  reaction  > 

0  +  f  CHr20  +  OH. 

The  calculations  were  insensitive  to  the  precise  value  used,  but  the  results  ar« 
consistent  with  recent  high  temperature  literature  values.  For  the  CHa,  mixture, 
the  best  agreement  was  achieved  using  a  recent  high  temperature  value  for 

0  +  Ctt*  t  CH-i  +  OH. 

The  data  indicated  that  this  reaction  has  a  markedly  non-Arrhenius  rate  constan' 
use  of  a  value  based  upon  low  temperature  data  gives  results  inconsistent  with 
those  observed.  The  CHi,  system  was  insensitive  to  any  other  methane  reaction 
rate  constant.  There  was  some  experimental  evidence  to  suggest  an  unsuspected 
complexity  in  methyl  radical  decay  channels  at  the  lowest  temperatures  observed. 

sThe  mechanism/ rate  constant  combination  used  here  was  then  applied  to  the 
analysis  of  literature  data  for  CH3  oxidation  by  O2.  This  analysis  was 
designed  to  obtain  a  value  for  the  rate  constant  of 

CH3  +  02  J  CH20  +  OH.  . 

Results  indicated  that  the  most  recent  high  temperature  studies  yield  values 
consistent  with  the  present  analysis. 
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